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Dynamical phase transitions in one-dimensional hard-particle systems 
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Department of Physics, University of Bath, Bath, BA2 7AY, United Kingdom 

We analyse a one-dimensional model of hard particles, within ensembles of trajectories that are 
conditioned (or biased) to atypical values of the time-averaged dynamical activity. We analyse two 
phenomena that are associated with these large deviations of the activity: phase separation (at 
low activity) and the formation of hyperuniform states (at high activity). We consider a version 
of the model which operates at constant volume, and a version at constant pressure. In these 
non-equilibrium systems, differences arise between the two ensembles, because of the extra freedom 
available to the constant-pressure system, which can change its total density. We discuss the re¬ 
lationships between different ensembles, mechanical equilibrium, and the probability cost of rare 
density fluctuations. 


I. INTRODUCTION 


Out-of-equilibrium systems differ from their equilib¬ 
rium counterparts in many ways. For example, long- 
ranged correlations are generic in non-equilibrium steady 
states [T1I2]; unusual phase transitions can take place [3]; 
and there may be significant differences in behavior for 
the same system in different ensembles (for example, 
canonical and grand canonical mm- Recently, there 
has been considerable interest in large-deviation phenom¬ 
ena [6], based on ensembles of trajectories that are condi¬ 
tioned on atypical values of time-averaged observables 13 - 
m- In glassy systems, ensembles conditioned on the dy¬ 
namical activity quite generically support phase transi¬ 
tions between active and inactive states mm]. In other 
systems, such biases can result in phase separation m, 
and other kinds of phase transition [3 . 

Here, we consider large deviations of the activity in 
a model of hard particles that diffuse in one dimension. 
Considering this system in a constant-volume ensemble, 
we showed previously m that it supports both phase- 
separated and hyperuniform states (in which large-scale 
density fluctuations are strongly suppressed M)- In con¬ 
trast to equilibrium systems, these biased ensembles of 
trajectories support coexistence of phases with different 
(mechanical) pressures: similar effects have recently been 
discussed in active matter systems [151 El]- We ana¬ 
lyze this effect using a Langevin equation for the time- 
evolution of the system - we explain it in terms of ran¬ 
dom forces that acquire non-zero averages in conditioned 
ensembles of trajectories. 

We also consider a version of the system in the 
constant-pressure ensemble. At equilibrium, the behav¬ 
ior in such a system resembles that of a subsystem of 
a large constant-volume system, and phase separation in 
constant-density systems is accompanised by bistable be¬ 
haviour in the constant-pressure ensemble. However, in 
the non-equilibrium ensembles considered here, this fa¬ 
miliar picture breaks down. We interpret this effect in 
terms of phase coexistence between states with different 
pressures. 

Taken together, our results highlight the broad range 
of phenomena that can occur in conditioned ensembles of 


trajectories, even in simple models. They show that equi¬ 
librium ideas of pressure and ensemble-equivalence can 
sometimes be applied in these contexts, but that these 
applications may be subtle and require careful justihca- 
tion. In the following. Section [H] describes the models 
and methods that we will use. Section |HI| shows results 
for constant-volume systems, while Section Jig includes 
results at constant-pressure. In Section |Vj we discuss 
force balance and mechanical equilibrium in these sys¬ 
tems, based around the virial pressure. We draw our 
conclusions in Section IVTl 


II. MODEL AND METHODS 
A. Model 


We consider a system of N hard particles undergoing 
Brownian motion in one dimension, with periodic bound¬ 
aries. Particle motion is described by the overdamped 
Langevin equation: 

= ( 1 ) 

kbT 

where Xi is the position of particle i, 17 is the potential 
energy of the system, and Dp is the diffusion coefficient 
for particle motion. The random noise rji has zero mean 
and is uncorrelated in time and space: 

( 2 ) 

We consider the case where the energy U is given by 
a pair potential, U = ^ ~ Specifically, we 

consider hard particles of size Iq, 


v{x) 


oo, a X < Iq 
0, otherwise. 


(3) 


In this case, particles cannot interpenetrate or other¬ 
wise move past one another. [For the derivatives in 0 
to make sense, one should regularize the potential by 
smoothing its discontinuities, so that the hard-particle 
case may be obtained by taking a suitable limit. In prac¬ 
tice, we will simulate the time evolution of Q using a 
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Monte Carlo scheme, so no explicit regularisation is re¬ 
quired.] 

In one dimension, the exact equation of state for this 
system is that of an ideal gas with excluded volume: 

P{L - Nlo) = NknT (4) 

where P is the pressure (we fix the units of energy by 
setting k^T = 1 so that the units of pressure are Iq^). 
We define the packing fraction (j) = 

We simulate this system in both the NVT ensemble 
(constant-density), and the NPT ensemble (constant- 
pressure). For simulations at constant-pressure P, the 
system size L evolves as 

L = -^[P+ (dU/dL)] + (t) (5) 

where is the diffusion coefficient for the volume 
coordinate. The noise term rj^ has zero mean, with 
{r]L{t)'i]L{t')) = S{t — t'), and ??L(t) is independent of the 
other noises r]i{t). 

To determine an appropriate value for Dl, we use a 
hydrodynamic argument. For large systems, the sys¬ 
tem size is almost always close to its mean value L = 
(L), with small fluctuations of size 0{\/N). Noting 
that {dU/dL) = “ ri)v'(ri+i - u) is the 

negative of the virial pressure, one may therefore lin¬ 
earise 0, taking [P + {dU/dL)] « ;^(T — L) where 

kt = {—l/L){dL/dP) is the isothermal compressibility. 
From the linearised Langevin equation, the volume re¬ 
laxation time is tll = {LnxkBT/D l). Diffusive scal¬ 
ing indicates that this relaxation time should be of order 
L /((27r)^Dc) where is a collective diffusion constant. 

Setting tll = L /((27r)^Dc) and assuming — Dp 
yields 

( 6 ) 

L 

The numerical prefactor in this equation is not crucial 
for this work but the scaling ol with system size L 
will be important in what follows. 


B. Monte Carlo dynamics 

We use a Monte Carlo (MC) method to implement 
these dynamical Langevin equations. The system evolves 
by single particle MC moves m- Particle displacements 
are chosen uniformly from the range —S < Ax < S', 
with S = O.l/o- If the displacement results in a particle 
overlap, it is rejected; otherwise, it is accepted. 

The natural unit of time in the system is tb, the time 
taken for the root mean squared displacement of a free 
particle to reach its own length: 


In terms of the Monte Carlo dynamics, this time interval 
corresponds to 


Nmc 


Nil 

S2/3 


( 8 ) 


attempted MC moves (here S^/3 is the mean squared 
displacement for a single step). This choice ensures that 
free-particle diffusion always has diffusion constant Dp, 
as required. The Monte Carlo dynamical scheme is equiv¬ 
alent to solving the Langevin equation Q, in the limit 
where the step size S' —>■ 0 [T5]. The value of S used here 
is small enough that qualitative features do not depend 
on S. 

When the packing fraction (f) is large, the rate of ac¬ 
ceptance of MC moves can get small, and it becomes 
convenient to use a rejection free MC algorithm that op¬ 
erates in continuous time [13 [19]. To achieve this, all 
possible particle displacements are calculated: let Qi de¬ 
note the fraction of possible moves for particle i that are 
compatible with the hard-particle interactions (i.e., the 
fraction that would be accepted). Particle i is selected 
with probability gi/9j) s-^id one of its possible moves 
is implemented. The simulation time is then incremented 

by 


ti = 


Tb 

Nmc 


^ -ln(l//i) 


(9) 


where /i is randomly distributed 0 < ^ < 1 (so ln(l//r) 
is exponentially distributed with a mean of unity). The 
process is repeated and moves are made until the total 
simulation time reaches the desired duration. 

In constant-pressure simulations, we also perform MC 
moves in which we propose changes to the volume of 
the system. A change in volume is proposed as Anew = 
Told + AL where —Sl< AL < Sl and the positions of 
particles are scaled uniformly by If this causes 

any particles to overlap the move is rejected immedi¬ 
ately. Otherwise, the move is accepted with probability 
min(l c 0 

At every MC step, a volume change is proposed with 
fixed probability 1/N and a particle displacement is pro¬ 
posed with probability {N — 1)/N. To satisfy equation 
([^ we specify that the maximum volume change satifis- 
fies = S‘^{2'K)^kBTKT/L. Note this depends on the 
applied pressure, via the mean box size L. 

We simulate high density systems where the typical 
free space per particle is much less than Iq and the mean 
collision time is much less than tb. In constant-density 
simulations we use (j) = 0.88; for constant-pressure we 
take the corresponding value P = 7.331^^. This gives 
a mean free space per particle of 0.136Zo and a typical 
collision time of 0.018rB. 


C. Large Deviations 


Tb = 


/2 

2D„ 


(7) 


As anticipated in the introduction, we will be con¬ 
cerned here with ensembles of trajectories which are bi- 
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FIG. 1: Equilibrium activity fc as a function of the packing 
fraction, (j), and the maximum step size, S, in the constant- 
density ensemble. The activity is linear in (f> for small packing 
fractions, because k is proportional to the number of particles 
per unit volume. For larger <j>, particles start to obstruct each 
other and the activity falls. 


tories of the system: 


Ps[x{t)] = 


iohs) 


( 12 ) 


where ^(s,tobs) = (e“'*^)eq is a dynamical partition 
function (the average is evaluated at equilbrium). By 
analogy with equilibrium thermodynamics, we interpret 
■0(3, fobs) = -{Ltohs)~^ log .^(s, fobs) as a dynamical free 
energy (or free energy density). For large N and tobsj 
the free energy may develop a singular dependence on s, 
which signals the presence of phase transitions. 

We use a notation (•).? to indicate averages with respect 


to the distribution (12). It is easily verified that the 
mean activity fc(s,tobs) = {Ltobs)~^{K)s is obtained by 
a derivative of the free energy: k{s) = {d'lp/ds), while 
the susceptibility x(s,tobs) = {Ltobs)~^ {{K - {K)s)'^)s is 
equal to —{dk/ds). 


ased to non-typical values of an ‘activity’ parameter K. 
Our analysis follows that of [9]. A trajectory of the sys¬ 
tem is a realisation of the system developing through 
time. A trajectory has a time duration fobs which is 
composed of M smaller and consecutive “slices” of time, 
each of length At such that fobs = At At. 

To define the activity K, we measure the squared dis¬ 
placements of all particles during each slice and sum over 
the M slices in a trajectory: 

M N 

( 10 ) 


D. Transition Path Sampling 


To sample the biased ensemble of trajectories given in 

ni. An initial 


(12), we use transition path sampling[] 


trajectory is generated from an equilibrium initial con¬ 
figuration by simulating the dynamics of a system for a 
duration fobs- A new trajectory is created by copying a 
randomly selected portion of the first trajectory to either 
the start or end of the new trajectory. The rest of the 
new trajectory is then generated according to the relevant 
Langevin equations. The new trajectory is compared to 
the first and replaces it with probability: 


-Paccept = min {l,e (13) 


where the notation [a:(t)] indicates a functional depen¬ 
dence on all particle positions throughout a trajectory. 
The quantity Axj = (l/N) - Xi ((j - l)At)] 

is the change of the centre of mass during the jth slice: 
for large systems this has a negligible effect on K but 
subtracting it in this way helps to minimize finite-size 
effects in simulations. It is often useful to normalise the 
activity as 


where AK = — Generation then continues us¬ 

ing the most recently accepted trajectory as the parent. 
Using this method means that, after many iterations, the 
algorithm samples trajectories according to the distribu¬ 
tion defined in ([12 ). 


III. RESULTS - CONSTANT-DENSITY 


U^obs 

(For the constant pressure systems, we replace L hy L 
in this equation.) We take the time interval At = tb, 
so that K depends on movement of particles on length 
scales comparable with their size, consistent with pre¬ 
vious studies 0. Fig. g shows the dependence of the 
activity k on the packing fraction (p and the MC step 
size S. 

The large deviation formalism uses a biasing field s to 
investigate trajectories where fobs is large but the activ¬ 
ity k deviates significantly from its typical value. To this 
end, we define a probability distribution over the trajec- 


In this section we present results for biased ensembles 
of this hard-particle system in the NVT (constant den¬ 
sity) ensemble. All results are for the case (p = 0.88, 
as in m- As discussed in that work, we expect many 
aspects of the system’s behaviour to be independent of 
p. 

A. Phase separation for s > 0 

We first consider the effect of a bias s > 0, which leads 
to a phase transition in this system[T2l [131122]- Phase 
transitions are signalled by singularities in the free en¬ 
ergy ip{s), which appear only in the limit when both the 
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(a) (b) 




FIG. 2: The average intensive activity for const ant-volume 
systems at 0 = 0.88. (a) Behavior of the activity for a system 
of = 100 particles, varying tabs- (b) The effect of changing 
the number of particles in the system at constant fobs = ZOtb . 


(a) 





FIG. 4: Trajectories of a constant-density system at A'’ = 60, 
(p = 0.88 and fobs = 20rB. Blue boxes represent particles. 
Time runs horizontally, position vertically. Applied bias of 
(A) s = -1, (B) s = 0, (C) s = 1.25. 


Fig. [3] shows a finite-size scaling analysis of the tran¬ 
sition. For fixed fobs = SOtb and increasing N, we find 
s* ^ N~^, as indicated by the collapse of k{s) and the 
scaling of the peak in the susceptibility y. It seems that 
the theoretical prediction that s* ~ N' -2 [Hill] may 
be observable only for larger-A^ and/or larger fobs: we 
discuss this possibility in Section [III C| 


FIG. 3: Scaling of the transition with system size, (a) k{s) 
collapses onto a single curve when scaled by N. (b) The 
peak in the dynamic susceptibility occurs at a const ant-value 
of s*N. Inset: The height of y* increases with increasing 
system size as the magnitude of fluctuations in the activity 
increases with N. 


observation time fobs and the system size N are very 
large. Fig. shows the average activity k{s) for differ¬ 
ent system sizes and observation times. In particular. 
Fig. Ha) shows the effect of increasing fobs at fixed sys¬ 
tem size N = 100, while Fig. [^b) shows dependence on 
system size N, all obtained for large fobs = 20 tb. As in 
glass-forming systems, one observes a crossover at some 
s = s*, from active to inactive dynamics; the value of s* 
is positive and depends on both N and fobs- 

We note that for s —>■ oo, the system must arrive at the 
state with minimal propensity for activity, which should 
be the fully phase-separated state, where the particles 
in the dense cluster are all touching each other. It is 
therefore clear that phase separation must occur at some 
field s*. In fact, this phase transition can be predicted 
and analysed in detail in the framework of fluctuating 
hydrodynamics [Tl|23l|24], which predicts that the system 
will phase separate whenever {d'^/dp'^){K)eq < 0, and 
that s* ~ A^“2 tends to zero as the system gets large [T^ 

mill]. 


B. Structure of the phase separated state 

Fig.i shows trajectories for the biased ensemble at 
several values of s, representing active, equilibrium and 
inactive trajectories. At s = 0 the system is an equi¬ 
librium fluid of hard particles: “bubbles” of local free 
volume are seen throughout the trajectory, distributed 
randomly in time and space. For s < 0 the system has 
slightly higher than equilibrium activity and exhibits no 
obvious structural change from equilibrium (see however 
Fig. I^below). In the inactive phase (s > s*, Fig|^) the 
system has fully phase separated for the whole trajectory. 
Temporal boundary effects can be seen at the beginning 
and end of the trajectory where the phase separation de¬ 
teriorates. 

To probe the structure of these systems it is convenient 
to make a change of co-ordinates. The system is one¬ 
dimensional and the particles are hard, so the ordering 
of the particle co-ordinates is fixed: there is no “overtak¬ 
ing” . We number particles so that their co-ordinates are 
in an increasing sequence, and define new co-ordinates 
Xj = Xj — jlo which are also ordered in the same way. 
The Xj are co-ordinates of point particles in a system of 
size L' = L — NIq. At equilibrium, the positions Xj are 
uncorrelated - they represent positions of ideal gas par¬ 
ticles. For example, if we define a Fourier-transformed 
density Spq = and calculate the structure fac- 
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(a) (b) 



FIG. 5: Structural measurements in the constant-density 
regime with point particles at ^ = 7.33 with N = 120 for 
equilibrium and inactive systems, (a) Measurements of the 
one-body density p(x) for homogeneous and phase-separated 
states, (b) distribution of separations for the equilibrium and 
inactive phases. (The dotted lines s = 0.125 are very close to 
the s = 0 results.) 

tor 

S{q) = ^{5pq5p-q) (14) 

then we find S{q) = N/L' = , independent of q. 

Now define di as the separation between particle i and 
its right neighbour, 

d^ = Xi+i - Xi. (15) 

In order to investigate the one-body density profile as¬ 
sociated with the phase-separated state illustrated in 
Fig. B c), it is necessary to fix an origin. We accom¬ 
plish this by finding the largest ‘gap’ di in any configu¬ 
ration. We choose a random point within that gap, and 
we place the origin the maximal possible distance, L' 12^ 
from that point. Thus the origin almost certainly lies 
within the dense phase. The density of point-particles is 
then p{X) = Ylj ~ where Xj is now measured 
with respect to this new origin. We average p{X) to ob¬ 
tain the one-body densities shown in Fig. [^a). At equi¬ 
librium, the density profile is uniform, as expected (up to 
weak boundary effects that arise because the origin was 
constrained to lie in the largest gap). As s is increased, 
the phase separation in Fig. appears as a non-uniform 
density profile: a large peak appears in p(x), which cor¬ 
responds to the clustering of particles near the origin. 

The distribution of separations, P{d), was also 
recorded for systems at s = 0 and s > s*. This is shown 
in Fig. Bb)- At equilibrium, one finds an exponential 
distribution, typical of a Id equilibrium fluid. However 
the distribution of separations in the inactive phase is 
bimodal. For small d, the distribution is approximately 
exponential but with a smaller characteristic length scale 
than the equilibrium fluid: this corresponds to particles 
within the dense region of the system. For larger d, there 
is a broad distribution of separations that comes from 
pairs of particles located on opposite sides of the large 
void (each configuration contributes only a single sample 


to the large-d peak, so the width of this peak appears only 
after averaging many configurations). These separations 
are comparable to the system size and are indicative of 
phase separation. 


C. Stability of the phase separated state 

The phenomenon of phase separation is not expected 
in one-dimensional systems at equilibrium (assuming 
that forces are short-ranged). We therefore explore the 
physics behind this effect more detail. One can estimate 
the probability that the trajectory shown in Fig. Qc) 
would occur at equilibrium, as follows. The particles 
within the large cluster are constrained by their neigh¬ 
bours and their contributions to the activity K are nec¬ 
essarily small. The particles at the boundaries of the 
cluster cannot move into the cluster, but they may move 
away from it. The probability that one of these particles 
nevertheless remains close to the edge of the cluster for 
the whole time fobs scales as where 7 is a pa¬ 

rameter with units of time, proportional to the rate for 
diffusion of the particle away from the cluster. 

The key point is that maintaining the integrity of the 
cluster requires only that the two particles at its bound¬ 
aries do not move away. Hence, for large fobs, the proba¬ 
bility Pps of a phase-separated state at equilibrium should 
satisfy ln(Pps/Peq) ^ ~ 27 tobs; where Peq is the prob¬ 
ability of a typical equilibrium trajectory. Within the 
s-ensemble, the ratio (Pps/Peq) is multiplied by a term 
e^^^, where AK is the difference in activity between 
equilibrium and phase separated states. Since the activ¬ 
ity at equilibrium is extensive (that is, AK ^ SkNtohs), 
one therefore expects phase-separated states to dominate 
for s > 2"f/{N5k). 

This argument essentially reproduces the prediction of 
[25] for phase transitions in kinetically constrained mod¬ 
els; see also m- Assuming that 7 is independent of N, 
we predict s* < N~^, consistent with Fig. As dis¬ 
cussed above, the more refined analysis available from 
fluctuating hydrodynamics predicts s* ~ N~'^ [T^]. That 
is, the bias required to stabilise phase separated states 
is even less than that predicted by the simple argument 
given here. The reason is that the interface between high- 
and low-density regions of the system may not consist of 
a single particle, but can be smoothed out: this acts to 
reduce 7 to a quantity of order 1/N, further stabilising 
the phase separated state. However, it seems that the 
regime in which these smoothed out interfaces can be 
observed is not accessible within our simulations (where 
both N and fobs are limited by the computational effort 
required.) 

D. Hyperuniformity for s < 0 

For s < 0, the system is biased so that the parti¬ 
cles move around more than they do at equilibrium. 
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(a) 



X 


found in jammed sphere packings [271I2S] and in a range 
of other physical systems [I31E9H3I]. 

In numerical simulations, the size of the system limits 
our ability to investigate the small-q behaviour of S(q): 
we can obtain data only down to a minimum wave vec¬ 
tor gmin = As shown in Fig. ^ density fluctuations 
at gniin are increasingly suppress^ as s becomes more 
negative. On increasing the system size, gmin is reduced, 
and S(qinin) becomes smaller (for any given s). Fig. |^b) 
shows that all values of S(q) collapse onto one curve at 
fixed s. This suggests that for any s < 0 then: 


(b) 



lim 5'((7niin,s) = 0 (16) 

N—>’OC 

while for s = 0 one has limAr_,.oo 5'((7niin, 0) = 

This discontinuous response to the field s at s = 0 corre¬ 
sponds to a dynamical phase transition to a hyperuniform 

state[T3]. 


IV. RESULTS - CONSTANT-PRESSURE 
A. Inactive state, s > 0 


FIG. 6: (a) The pair correlation function, g{x), for point par¬ 
ticles in the constant-density regime with N = 120, fobs = 
20rB in the active phase. The depression in g{x) for small x 
reflects an increase in local free volume around the particles 
as activity increases, (b) Small q structure factor measure¬ 
ments for constant-density systems. Biasing to s < 0 causes 
the suppression S{q) ~ g consistent with the onset of hyper¬ 
uniformity. At fixed s all system sizes collapse onto a single 
curve. 


Fig. ^a) shows the pair correlation function g{x) = 
{p{x')p{x'+x))/{p)'^ in this regime. (The system is trans- 
lationally invariant so there is no dependence on a:'). Par¬ 
ticles appear to repel each other, leading to a depletion 
zone around each particle, which facilitates motion on 
small length scales. However, the deviation of g{x) from 
unity is small in absolute terms: this is a rather weak 
effect. 


On the other hand, moving to Fourier space reveals a 
much stronger response to the bias s, associated with 
long-ranged correlations in the density. This is illus¬ 
trated in Fig. [^b), which shows the structure factor 
of the system for s < 0. At small wavevectors (large 
length scales) the structure factor shows markedly dif¬ 
ferent behaviour from the equilibrium fluid. At equi¬ 
librium the system behaves like an ordinary fluid with 
S{q) = const, for all q [26], recall (14|. Trajectories for 
s < 0 have a very different structure factor as q ^ 0. 
This suppression of long range density fluctuations is a 
sign of hyperuniformity [14]: Hyperuniform states are dis¬ 
tinguished by anomalously small density fluctuations on 
long length scales, eventually vanishing at infinite range, 
hence S{q) —)■ 0 as g —>■ 0. The phenomenon has been 


We now consider the constant-pressure version of the 
model, in which the system size evolves in time according 
to equation ([^. In equilibrium, one expects local prop¬ 
erties of single phases to be independent of ensemble. 
For example, since the pressure is constant throughout 
an equilibrium system, one may think of the constant- 
pressure simulation as representing a subsystem of a very 
large constant-volume system. One might expect the 
same equivalence to hold in biased ensembles at s 0, 
but we will see that the applied bias s leads to a break¬ 
down of ensemble-independence. 

We take the pressure P = 7.331q^ so that the mean 
volume fraction at equilibrium is (0) = 0.88, consistent 
with Sec. jlllj The effect of the biasing field s is shown in 
Fig. [7] Comparing with Fig. [^ a similar transition is ap¬ 
parent, but instead of a crossover at s* > 0 that depends 
on N,tohs, one instead observes a crossover very close 
to the equilibrium point s = 0. Note that the values of 
fobs used here are significantly larger than those used in 

the constant-volume system: they are comparable with 

_2 

the volume relaxation time of the barostat ^ L (re¬ 
call Section HB |. Fig.j^shows representative trajectories 
from biased ensembles in the constant-pressure system. 
Comparing with Fig. [^ no phase separation occurs. We 
also note that the box size varies with s, consistent with 
Fig. 00). 

Fig. d shows the correlation between activity and 
global density for all dynamic regimes. Within fluctuat¬ 
ing hydrodynamics, the density field is assumed to give 
a full description of the large-scale behaviour of this sys¬ 
tem, In constant-volume systems, density fluctuations on 
finite wave vectors control the fluctuations in activity 
[Hdl]. In the constant pressure system. Fig. shows 
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(a) (b) 



(c) 



FIG. 7: Dynamic behaviour of the constant-pressure system. 
All results come from trajectories of duration tabs ~ at 
P — 7.33. (a) the intensive activity and dynamic susceptibil¬ 
ity of the system as a function of s. The transition (crossover) 
takes place at s* = 0 for all system sizes, and the width of 
the crossover is proportional to (b) The data for all 

system sizes collapses when scaled by tll, (inset) the peak 
in the dynamic susceptibility scales with tll- (c) The (total) 
density of the system changes as the system undergoes the 
transition shown in (a). 



Time 


FIG. 8: Representative trajectories of constant pressure sys¬ 
tems at P = with N — 40, fobs = 220tb at different 

biases. Red boxes represent the boundaries of the system. (A) 
s = —0.250, the active systems increase in length and thus re¬ 
duce global density. (B) s = 0, the equilibrium system has a 
fluctuating volume but maintains {(p) = 0.88. (C) s = 0.375, 
the inactive system is compressed relative to equilibrium and 
thus has suppressed activity. 
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FIG. 9: Scatter plot of activity against the average density of 
trajectories of duration 1300rB with N = 100 particles from 
biased ensembles. There is a strong correlation between the 
activity and the density of the trajectories. 


that the total density correlates strongly with the activ¬ 
ity: high density is associated with low activity, and vice 
versa. 

Given this strong correlation, we can link the phase 
transition that takes place at s = 0 in this system to 
the diverging hydrodynamic time scale tll (recall Sec¬ 
tion IIA). We define the (normalised) autocorrelation 


function of the system size 


ChLit) = 


{SL{t')SL{t' + t))eq 


mm 


(17) 


eq 


which is evaluated at equilibrium (so there is no depen¬ 
dence on t'), with SL = L — L. This correlation function 
decays on a time scale close to tll- Similarly the corre¬ 
lation function of the activity is 


(<5/C(P)2)eq 


(18) 


where lC{t) = \ri{t +l2^t)—ri{t) — Aa;(<)p is the quan¬ 
tity that appears in the definition of the activity K, re¬ 
call ( |Io| ). To show the long-time behaviour of Ckk{t) 
more clearly we smooth the function by convolving it 
with a Gaussian window, with variance cr^ = Tg/4: we 
plot where the proportionality 

constant T normalises the correlation function to unity 
at t = 0. 

The correlation functions Ckk{t) and CkLit) behave 
very similarly, consistent with the idea that the activ¬ 
ity fluctuations are strongly correlated with those of the 

global density (and hence to the system size). Since the 

_2 

volume relaxation time tll diverges as L , we therefore 
expect a similar divergence in the relaxation time of the 
activity. 

This divergent time scale is important because the sus¬ 
ceptibility X is related to the autocorrelation function of 
the activity as 

2 

x{s = 0, fobs ^ oo) = = / dt ((5/C(0)(5/C(t))eq (19) 

^ Jo 
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FIG. 10: The activity and volume correlation functions for 
constant pressure systems at equilibrium, for different system 
sizes. The correlation time, Tkk, is strongly correlated with 
tll- To display the long-time behavior most clearly, Ckk{t) 
has been smoothed with a Gaussian window (see main text). 


_2 

so that X ~ '^kk ~ L diverges at s = 0, which we inter¬ 
pret as a dynamical phase transition. (The equal time 
value of the correlator in this equation scales as L since 
61C is extensive in the system size: this L-dependence 
cancels with the prefactor so that the right hand side 
scales with Tkk, with a prefactor of order unity.) Since 
y = —dk{s)/ds, the divergence of x corresponds to a 
singularity in fc(s) and hence a dynamical phase transi¬ 
tion. This amounts to a perturbative argument for the 
existence of the phase transition: a related perturbative 
argument based on fluctuations at finite wavevector was 
used in [13| to explain the existence of phase transitions 
in systems at finite volume. 

Recalling the discussion of Sec. |III C] the analogous ar¬ 
gument for the constant-pressure system is that the noise 
force in ([^ acquires a finite average in the inactive 
state, resulting in a reduced system size. The previous 
argument based on indicates that biasing the noise 
force in this way requires very little cost in probability: 
this low cost appears partly because only a single noise 
term needs to be biased, but also because the absolute 
size of the bias becomes small in large systems, due to 
the scaling of the diffusion constant Dl with system size 
(Eq. (§). 

Finally we note that s* shifts from a value of order 1/A^ 
in the constant-density system to a value close to zero at 


(b) 



FIG. 11: (a) Structure factors in the constant-pressure regime 
when biased to higher than equilibrium activities. There is 
a suppression of long range density fluctuations as in the 
constant-density active phase, (b) The distribution of sep¬ 
arations for a system of = 100 particles. Symbols repre¬ 
sent measured distributions, solid lines are exponential dis¬ 
tributions with a mean separation calculated from the mean 
volume. In all regimes the separations are distributed expo¬ 
nentially and are similar to equilibrium, albeit with a different 
mean separation. 


constant pressure. We interpret the small positive s* (for 
the const ant-volume system) in terms of the probability 
cost required to form the interface in a phase-separated 
system at constant density. At constant pressure, no in¬ 
terface is required so the system can have a diverging 
linear reponse, as shown in (19). 


B. Active state, s < 0 

In Fig. [TT|(a), we show the structure factor of the 
constant-pressure system for s < 0. For a given bias 
s, the fluctuations in the total system size are small in 
relative terms, so we evaluate the structure factor at 
wavevectors q = 2m:jV as usual, and calculate S{iq) by 
an ensemble average at fixed n. This provides an esti¬ 
mate of S'(( 7 ) with q = 21 : 11 /L'. The results of Fig. [ll|(a) 
are consistent with hyperuniformity of the active (s < 0) 
phase, although the effect is weaker than that shown in 
Fig.i for the constant-density system. We also show 
the distribution of particle separations in Fig. Sb) , for 
comparison with Fig. [^b). The distribution fits well to 
an exponential form, independent of s. Given the corre- 
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lations that are apparent from Fig. 11 |a), this result is 
somewhat surprising: it might be that the correlations 
are sufficiently weak on short length scales that they are 
not discernible from P{d). 

In this constant-pressure system, it seems that achiev¬ 
ing large deviations by changes in structure (for exam¬ 
ple phase separation or hyperuniformity) is unfavourable 
compared to changing the system density. Thus, par¬ 
ticle separations remain exponentially distributed inde¬ 
pendent of s, but the system density depends strongly on 
s. Even at the longest length scales and largest |s| the 
structure is only weakly affected, compare Fig. [^a) and 
Fig.^a). 


V. DISCUSSION - PRESSURE BALANCE AND 
MECHANICAL EQUILIBRIUM 


distance Lb then = —LbPb^'^'^: where the overbar 
denotes a time average (over fast degrees of freedom) and 
pmech jg average magnitude of which is the me¬ 
chanical part of the pressure applied to the region B by 
its environment. 

In addition, recalling Vb = the equation 

of motion 0 yields 

Vb = ^ x.idtXi - r]i) 

= ^ [5t(a;?/2) - x,r]i] (23) 

P iinB 


For a system in a steady state, one has J2iinB ~ 

and assuming that fji = 0 then Vb ~ 0. In this case ( 22 ) 
yields 


To rationalise the numerical results shown here, it is 
useful to consider the balance of forces in the phase- 
separated state shown in Fig 0 - To this end, we con¬ 
sider the virial formula for the pressure. As we briefly 
review here, this expression can provide information on 
force balance even in systems that are not at equilibrium. 
We define the virial by 

V = ^a:,F, (20) 

i 


where Fi is the conservative (non-stochastic) force on 
particle i, which in our system is Fi = v'{xi+i — Xi) + 
v'(xi-i — Xi), where v(x) is the interparticle potential (re¬ 
call ([^) and the prime denotes a derivative. (An alterna¬ 
tive would be to include both conservative and stochastic 
forces in the definition of V, but we choose this definition 
for later convenience.) 

By restricting the sum in (20) to particles within a 
given region of the system, and time-averaging over a 
period long enough for fast degrees of freedom to relax, 
we now show that we can obtain an estimate of the local 
pressure. Our discussion follows that of [32] ■ We con¬ 
sider a region B of the system, of linear size Lb- For 
convenience, we assume that this region does not cross 
the periodic boundaries of the system. Then, if Vb is the 
contribution to (pOl) from particles i in B then 


P imech 
B 


v^ 

Lb 


(24) 


which is equivalent to the standard result for the mechan¬ 
ical part of the virial pressure. (In the usual case [32] , one 
considers Newtonian dynamics instead of overdamped 
Langevin dynamics, in which the pressure includes an 
extra idea gas term Nk^T/L b- For overdamped dynam¬ 
ics, this term does not appear, but this has no impact 
on mechanical equilibrium because the ideal gas contri¬ 
bution to the pressure is exactly constant throughout the 
system.) The approximate equality in (24) appears be¬ 


cause we have invoked a time-average over fast degrees of 
freedom, instead of a full ensemble average. This distinc¬ 
tion is useful when applying the result to phase separated 
steady states such as that shown in Fig. 0 . 

To consider those states, the preceding argument must 
be generalised, because the biasing field s can lead to 
noise forces rji with non-zero averages, iy, 7 ^ 0. In 
this case, balance of the total forces on region B yields 
Apmech ^ ^ where AP“®®1^ is the difference 

in mechanical pressure between the two sides of the re¬ 
gion. Considering a small region B centred at position 
cc, we obtain a force-balance equation 


VP’ 


mech 


(x) = p'q{x) 


kBT 

L>n 


(25) 


Vb = ^ + 6j^i_i)v'{xj - x,)- (21) 

2 in S j 

If we also decompose the sum over j into contributions 
from inside and outside B then we obtain 

Vb = Vt‘ + V|r‘ (22) 

with = ^Y.^jinBixj - Xi)v'(xi - Xj)Sj,^+i and 

V°B^ = SiinB where is the force on parti¬ 

cle i from those particles outside B- In a steady state, 
the forces on the two boundaries of the region B must 
balance. Since the two boundaries are separated by a 


where the noise force is averaged over particles in the 
vicinity of position x- Again for small boxes (of size 

5x), one finds V{x) = —’^^xSxpfjix) and V'^’'*(x) = 
-feP“®®h(x) - a:feVP“®®‘'(x). Hence from ( |22|25| ): 


P’ 


mech 


(x) = 


i;int( 3 .) 


5x 


(26) 


That is, the virial still gives a useful estimate of the local 
pressure, even in the presence of noise forces with non¬ 
zero averages. Note however that in the presence of these 
noise forces, the mechanical pressure is not constant in 


space, but varies according to (25). 
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Returning to the phase separated state in Fig. [^, we 
see that the dense and sparse regions of the system are 
both homogeneous, so we expect the local virial V“*(a;) 
and the mechanical pressure to be constant within each 
phase. However, it is easily verified that the virial dif¬ 
fers strongly between the phases, so the mechanical pres¬ 
sure also differs. The origin of this pressure difference is 
the presence of non-zero noise forces at the boundaries 
of the cluster, consistent with (251 and the discussion 
of Sec. Ill C Given this observation, it is not surprising 
that phase coexistence was not observed in the constant 
pressure system for s > 0: the phases that coexist in the 
contant-volume system do not have equal pressures. 

We note in passing that this argument may be straight¬ 
forwardly generalised to active matter systems [El [16], in 
which case the noise forces rji should be replaced by the 
particles’ self-propulsion forces. In scalar active matter, 
one may find coexistence between states with different 
mechanical pressures, due to the presence of one-body 
forces with non-zero averages. 


VI. CONCLUSION 

We have demonstrated dynamic phase transitions in 
a one-dimensional model of diffusing particles, including 
transitions from simple equilibrium fluid states into both 
high-activity and low-activity states. We considered both 
constant-density and constant-pressure systems: their 
transitions share some common features but there are 
also important differences. Based on the theory of fluc¬ 
tuating hydrodynamics, we argued m that the transi¬ 
tions occur for all densities p, and the arguments given 


here indicate that this should also hold for all applied 
pressures. 

Considering first the transition to inactive states, 
the constant-volume system undergoes phase separa¬ 
tion, while the constant-pressure system increases the lo¬ 
cal density. (As in equilibrium systems, the constant- 
pressure system avoids interfaces between coexisting 
phases.) For large systems, the inactive phase occurs 
for all s > 0, in both ensembles. However, there is no 
signature of phase coexistence in the constant-pressure 
system: we find only a dense phase, consistent with the 
different (mechanical) pressures of the coexisting phases 
in the contant-volume system. 

For transitions to high-activity states, the constant- 
volume system spontaneously suppresses long range den¬ 
sity fluctuations and develops a hyperuniform struc¬ 
ture m- At constant pressure, the main feature of 
the high-activity state is that the total density decreases 
sharply, although there is also some suppression of large- 
scale density fluctuations. 

Overall, these results emphasize that equilibrium ideas 
of ensemble equivalence do not apply directly when con¬ 
sidering large deviation phenomena such as those consid¬ 
ered here. While the mechanical pressure and the virial 
can still be related, the possibility of phase coexistence 
at unequal pressures shows how familiar equilibrium con¬ 
cepts such as phase separation need to be re-evaluated 
and generalised in these non-equilibrium settings. 
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